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Solutions for the static spherically symmetric extremally charged dust in the Majumdar- 
Papapetrou system have been found. For a certain amount of the allocated mass/charge, the 
solutions have singularities of a type which could render them physically unacceptable, since the 
corresponding physically relevant quantities are singular as well. These solutions, with a number 
of zero-nodes in the metric tensor, are regularized through the (5-shell formalism, thus redefining 
the mass/charge distributions. The bifurcating behaviour of regular solutions found before is no 
longer present in these singular solutions, but quantized-like behaviour in the total mass is ob¬ 
served. Spectrum of regularized solutions restores the equality of the Tolman-Whittaker and ADM 
mass, as well the equality of the net charge and ADM mass, which is the distinctive feature of 
Majumdar-Papapetrou systems. 

PACS numbers: 04.40.Nr 


I. INTRODUCTION 

In this paper we present a class of newly found singular, static, spherically symmetric solutions to the Einstein- 
Maxwell equations for the extremally charged dust (ECD) model of matter. The solutions for the metric profile 
functions show remarkable discrete properties with respect to the amount of allocated mass. These metric profile 
functions have nodes that will be shown to correspond to naked singularities and they are mimicking the behaviour of 
a quantum system with regard to its higher and higher excitations. Regularization procedure for the removal of the 
naked singularities that are present in our solutions is investigated. The regularization is accomplished by introducing 
a thin shell (<5-shell) of ECD into the system. The regularized solutions retain their properties with respect to the 
amount of the allocated mass, as well as their asymptotic behaviour (charge and ADM mass). 

We first briefly review the equations of structure in the Majumdar-Papapetrou (MP) formalism [1-3]. In the next 
section we present the new solutions with non-negative energy density of ECD. These (singular) solutions exhibit zero 
crossings (nodes) in the metric profile function. In Sec. Ill we show that introduction of a spherical (5-shell of ECD 
suffices to remove the spacetime singularity. Conclusions are given in Sec. IV. 

The system is determined by the coupled Einstein-Maxwell equations, 

R i1 v -^R = S,'kT^ v and = Anf, (1) 

where is the energy-momentum tensor, F^ is the electromagnetic field strength tensor, and j /J is the electric 
charge 4-current. Denoting the 4-velocity of the source and a the charge density, we have j /l = au^. The 
energy-momentum tensor due to the electrically charged (pressureless) dust, T^ v , consists of the matter component 

( 2 ) 

where p is the energy density of the dust, and the electromagnetic component 

E ^ = ^(- F ^ F ™ + l^F^F pa ) (3) 
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due to field strength F^ v = — d v A tl . As it was shown by Majumdar [1] and by Papapetrou [2], if the energy 

density of the static dust equals the charge density, p = ±<r, one can write the line element using the harmonic 
coordinates, 

ds 2 = U~ 2 dt 2 -U 2 (dX 2 + dY 2 + dZ 2 ), (4) 

where U = U(X, Y, Z ), and the Einstein-Maxwell equations (1) for the energy momentum tensor + E^ v 

reduce to a nonlinear version of the Poisson equation, 


V 2 f/ = -4irpU 3 (5) 

(V 2 is the 3-dimensional flat space Laplacian operator). Such matter is called the extremally charged, or electrically 
counterpoised dust (ECD). The 4-velocity of the static ECD is = U Sq , the charge 4-current is = aU Sq, and 
the 4-potential has only the time component, A M = d°<f>, where the electrostatic potential $ is related to the metric 
profile function by $ = ±([/ -1 — 1). For a formulation of the MP formalism in d > 4 dimensional spacetimes, as well 
as for a recent compilation of references dealing with MP systems, see [3]. 

While one of the most remarkable features of the MP systems is that they do not require any spatial symmetry, 
spherically symmetric and asymptotically flat MP systems have played a role in obtaining valuable insights into the 
properties of astrophysically plausible systems, and are sometimes called ‘Bonnor stars’ [4], In the spherical symmetry, 
the harmonic line element (4) is 


ds 2 = U(R)~ 2 dt 2 - U(R) 2 (d R 2 + R 2 dfi 2 ), (6) 

where df2 = sin 2 d dip 2 + d$ 2 , and the MP equation (5) can be written 

R~ 2 (R 2 U'Y = -dirpU 3 , ( 7 ) 

the prime denoting differentiation with respect to the harmonic coordinate R. The energy momentum tensor due to 
the electrostatic field (3) acquires the diagonal form: 

= <f>' 2 diag(+,+,-,-)/87 t. (8) 

The total energy density, p t ot. = T °o = p + <E> ,2 /87t, involves the contribution from the ECD which is non-negative 
by assumption, and the manifestly non-negative contribution from the electrostatic field. The pressures are of purely 
electrostatic origin. The transverse pressures, p tI . = —= 4> ,2 /87t, are non-negative, while the radial 
pressure, p ra( j. = —T r r = — $ ,2 /87r, is of equal strength, but of the opposite sign (tension). These relations can be 
summarized as follows: 


Ptr. — Prad. d 5s Ptr. ^ Ptot. — Ptr. T P- (9) 

The weak, null, strong, and the dominant energy conditions [5] are all satisfied, as expected of classical matter and 
fields. 

The Maxwell equation can be written in terms of the Hodge-dual of the field strength tensor, the 2-form *F a p = 
and the Hodge-dual of the charge 4-current, the 3-form *j ot p 1 = , as d*F = dir*j. By the 

fundamental theorem of exterior calculus (Stokes’s theorem), one can use the Maxwell equation to express the charge 
Q contained within a compact 3-region 1Z as the integral of the field over the 2-boundary dlZ, Q = J lz *j = f n d*F = 

A- f dn *F. In our context, we use the t = const, slice as the region 1Z, and the charge can be written as 

Q = (io) 

where R = R\ and R = R2 are the radii of the inner and the outer spherical 2-surfaces acting as boundaries of the 
3-region 7 Z (we have used j° = <rU and cr = ±p = =FV 2 I//47rt/ 3 ). In the regular center the electrical field vanishes 
(i.e. = U' = 0 at R = 0), and as i?i —> 0, the flux of the field through the inner boundary vanishes as well. In 
fact, the inner boundary can be omitted, and the total charge of the system is obtained by letting Ri —> 00 , giving 

Q = T hniR^oo R 2 U'. 

The MP geometries are static and therefore they allow for the timelike Killing vector field satisfying = 0, 

orthogonal to the t = const, slice. Since the geometries are also asymptotically flat, one normalizes the Killing 
vector so that liiiiR^Qc = 1- The mass M can be given by the Komar’s formula [6] as the integral of the 
2-form a a p = — e Q / 3 A t ! yV At £ 1 ' over the 2-sphere S of sufficiently large radius, M = f s a, which coincides with the 
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ADM mass of the system. However, by the Stokes’s theorem, this integral can be expressed as the integral of 
the 3-form da (the exterior derivative of a) over the interior of the 2-sphere, f int ^ da. It can be shown [7] that 

(da) a/ 3 7 = — 47 e a/ 3 7P V"V ! ,£^ = j^e a / 3 1 ij,R fJ ’ v £, v = e a /37^(2 T^ v — . In the case of the static spacetimes, this 

reduces to the Tolman-Whittaker expression [8, 9], which in our context integrates to 


M = J (2T° 0 - T)y/=g d 3 A 


U' 2 
47rC/ 4 


V 2 [/\ 

AnU 3 ) 


V^d 3 X 


R 2 U' r 2 
U R 1 ’ 


( 11 ) 


where R\ and i ?2 are the radii of the inner and the outer spherical surfaces (boundaries). In a geometry with the 
regular center we can again omit the inner boundary, and the mass of the system is obtained by letting i ?2 —> oo, i.e. 
M = - liniR^oo R 2 U'/U. 

Therefore, in asymptotically flat and nonsingular MP spacetimes we have 


M = ±Q = (ADM mass), (12) 

where the mass M and the charge Q can be given by the integrals over the t = const, slice, and their values coincide 
with the value of the ADM mass obtained at spatial infinity. However, when singularities in the metric are present, 
the ‘volume’ integrals (10) and (11) cannot be straightforwardly applied. We deal with such solutions in the next 
Sections. 

It is also important to give the relation among the harmonic coordinates (6) and the curvature coordinates with 
the line element 


ds 2 = B(r) dt 2 - A{r) dr 2 - r 2 dH 2 , (13) 

where r = UR, and the metric components are related by A -1 / 2 = 1 + ( dU/dR ) R/U, and B = U~ 2 . 

In absence of ECD, the general spherically symmetric solution of (5) is U{R) = k + m/R 1 where k and m are 
integration constants. To obtain asymptotically flat spacetime at infinity one sets k = 1, so in the spherically 
symmetric harmonic coordinates (6) the metric profile function is 

U(R) = l + m/R, (14) 

while in the curvature coordinates (13) the metric profile functions read 

B{r) = A -1 (r) = (1 — m/r) 2 . (15) 

For m > 0 this is the extremal Reissner-Nordstrom (ERN) black hole spacetime of ADM mass m. At the harmonic 
radius R = 0 (corresponding to the curvature radius r = m), there is the event horizon, while at R = —m (r = 0), 
there is a point-like singularity. In the case of negative ADM mass, m < 0, at R = \m\ (corresponding to r = 0) there 
is a point-like singularity, and since there is no event horizon in the spacetime, this singularity is naked. 


II. THE SOLUTIONS 


The MP solutions can be constructed for the assumed distribution of the ECD specified by the non-negative function 
p(R ). If p(R) vanishes as R —» oo more rapidly than R~ 3 , the solutions are asymptotically flat at infinity. This is 
expressed by the boundary conditions 

lim U(R) = 1 and lim R 2 U'{R) = —moo ■> (16) 

R—>oo R —xx) 

where jtt-oo is the ADM mass of the system. 

We will be looking for the solutions in which the metric profile function U is regular (finite) at R = 0. In the 
solutions in which U has no nodes, i.e. it is finite for all R > 0, the resulting spacetime will be free of singularities. If, 
on the other hand, U(R) has one or more nodes at R > 0, then the outermost node corresponds to a naked pointlike 
singularity. That the spacetime is singular can be seen by observing that the Riemann quadratic invariant, 


R^ kS R^ kS = |l4 R 4 U' 4 + 2 R 2 U U '- 2 (2 RU' — 5 R 2 U") 


+ U 2 (8R 2 U' 2 +4R 3 U'U" + 3R 4 U" 2 ) , (17) 
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evidently diverges if for some Rq > 0, U(Rq) = 0. The Ricci curvature scalar, 

R= (4 U' + 2 RU")/(RU 3 ), 


(18) 


and the Weyl quadratic invariant, 


C^ kS C^ kS = 16 [3RU' 2 + U (JJ' — RU")] 2 /(3R 2 U 8 ), (19) 

also diverge as R —> Rq. The singularity is pointlike since the curvature radius r, which is related to the radial 
coordinate R by r = RU(R), tends to zero as R —» Ro- The singularity is naked because there is no event horizon to 
protect it from the outer space. 

We will model the distribution of ECD with the exponential function 

p{R) = VPo exp [-(R/R) 2 }, ( 20 ) 

where 77 is the adjustable source strength factor, and po and R are constants, i.e. we will solve 

R- 2 {R 2 U')' = -4irr]p(R)U 3 (R). (21) 

It was found before [10] that the regular solutions do not exist for values of 77 greater than some critical value rj c , while 
in the range 0 < 77 < r] c , the critical solution bifurcates into two independent solutions with different ADM masses. 
The ‘high mass bifurcation branch’ produces spacetimes that asymptotically coincide with the external part of the 
ERN metric, they allow arbitrarily large red-shifts [4, 11], and are known as ‘quasi black holes’ [10, 12]. 

We obtained the solutions to (21) for the metric profile functions U(R) with zero, one or more nodes, subject to 
the boundary condition (16) with the fixed value of (ADM mass). At R = 0, we required that the metric profile 
function is finite (see the above analysis), and we found that such solutions exist only for a discrete set of values of 
the source strength factor ry, i.e. for discrete amounts of mass allocated in the source term of (21). The ‘quantized’ 
behaviour of the metric profile functions with nodes is shown in Figs 1 and 2. 

Starting from the solution Uq(R) which is finite at R = 0, 0 < t/o( 0 ) < oo, and has no nodes, corresponding to 
source strength ryo, a slight increase/decrease in the source strength, ryo —* ryo ± Sr), makes U(R) diverge towards =Foo 
at R = 0 (see Fig. 1 upper plot). If we increase the source strength further, the divergence persists until a particular 
value of the source strength rji is reached, for which the solution Ui(R) is again finite at R = 0, but now with a 
negative value, —oo < C/i(0) < 0, thus U\(R) having one node (see Fig. 1 lower plot). In this way one can generate 
solutions with two, three, and more nodes (see Fig. 2). 

Since all these solutions obey the same boundary condition as R —> oo, they cannot be distinguished by a distant 
observer, i.e. they all have the same ADM mass moo- For the regular solutions (no nodes), the charge integral (10) and 
Tolman-Whittaker mass integral (11), integrated over the t = const, slice, obey the relation (12), M = ±Q = m^. 
In the case of the solutions with one or more nodes in the metric profile function, the outermost being at R = Rq, 
the integration of (10) and (11) over the t = const, slice from R = R 0 to R . = oo does not give the expected value 
moo- This can be seen by looking at the nonvanishing contributions of the fluxes through the infinitesimal inner 
boundary (the 2-surface of radius R —> Ro, corresponding to r —> 0) to the respective integrals. In the case of the 
charge integral (10), this term contributes the finite value Q o = ±R 2 U'\r 0 , thus rendering Q = ±moo + Q o- This can 
be understood as — Qo being the point charge located at the singularity. In the mass integral (11), the flux through 
the inner boundary diverges, leading to M = oo. In analogy with classical electrodynamics, this can be understood 
as the diverging electrostatic energy of the point charge. 


III. REGULARIZATION WITH THE .5-SHELL 

In view of the above results, we now look for the possible modification of the ECD distribution that can remove 
the pointlike naked singularity located at R = Rq (corresponding to the outermost node of the metric profile function 
U) while keeping the geometry outside of some radius R = a > Rq unchanged. In the same time we require that in 
the regularized solution the charge integral (10) and the Tolman-Whittaker mass integral (11), when integrated over 
the t = const, slice, obey the relation ( 12 ). 

All the singular configurations considered so far can be regularized by introducing a single spherical shell of ECD 
with the regular geometry in the interior. We will be using the (5-shells with flat interior geometry, expecting similar 
behaviour also for the shells of finite thickness [13, 14]. The flat geometry in the interior (with constant metric 
profile function U) is to be joined with the exterior geometry of our solutions, given by the metric profile function 
U(R) satisfying (21) and with the outermost node at R = Rq, along a spherical (timelike) hypersurface given by 
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U(R) 



FIG. 1: The metric profile function Uo(R), finite at R = 0 with no nodes (thick line, upper plot), and the function Ui(R) finite 
at R = 0 but with one node (thick line, lower plot). The metric profile functions with the same asymptotic mass, corresponding 
to 20% higher and 20% lower source strength r] and diverging as R —» 0 are shown for both cases (thin lines). Metric profile 
function for the ERN spacetime is shown for comparison (dashed line). 

U(R) 



FIG. 2: Metric profile functions Ui(R) with i = 0,..., 3 nodes, finite at R = 0, are shown. (Non-linear scaling is applied on the 
R axis.) 
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R = a > R 0 . Applying the Israel’s junction condition formalism [5, 15], we will obtain the components of the required 
energy-momentum distribution on the junction hypersurface (shell). On the both sides of the hypersurface we use 
the coordinates = (t, R, •&, if), while as the coordinates on the hypersurface we naturally adopt y a = ( t , $, if). 

The first junction condition requires that the induced metric on the hypersurface (the first fundamental form of 
the metric on the hypersurface), given by h a b = where e(] = dx^/dy a , must be the same as induced by 

the metrics on either side of the hypersurface, i.e. the discontinuity of the induced metric over the hypersurface must 
vanish. In our context, this requires that the metric profile function U is continuous over the hypersurface. The 
regularized metric profile function C/ reg can therefore be written 

U iee (R) = U(a) + (U(R) - /7(a)) 6 (R - a), (22) 

where 0(R) is the unit step function. 

The second junction condition requires that the extrinsic curvature (the second fundamental form of the metric on 
the hypersurface), given by K a b = V„n M , where n M = 5^U is the unit normal vector to the hypersurface, must 
be the same as computed from the metric on either side of the hypersurface. Here, the nonzero components of the 
extrinsic curvature are 


K° o = U'/U 2 , K** = K v v = -( U/a + U’)/U 2 


(23) 


where U and U' are evaluated at R = a. Since in the flat interior we have U' = 0, while in the exterior we have 
U’ ^ 0 in general, the second junction condition cannot be satisfied. However, given that the first junction condition 
is satisfied, it is still possible to join the two geometries, provided that we introduce (allow for) the appropriate 
distribution of the energy-momentum on the hypersurface (shell). The components of the energy-momentum tensor 
on the shell are then given by the expression 

^=^([[A' a6 ]]-[[A']]/i ab ), (24) 

where K = K° a and the double brackets denote the discontinuity of a quantity over the hypersurface. In our context, 
for the surface energy density of the shell we obtain 

a = S° o = —U’/A-kU 2 , (25) 

while the surface tension, 6 = S®# = S vanishes identically. From this one can see that the energy-momentum 
distribution on the shell satisfies all standard energy conditions if a is non-negative, while all energy conditions are 
violated if a is negative. Since in our solutions the ADM mass is positive, the metric profile function U is outwardly 
decreasing (16) as R —> oo, but since it is continous and nonnegative for R > Rq it must have a maximum at some 
R = 7/max > Rq- It follows that if the regularizing shell is introduced at R = a € {Ro,Rmax), he. where U is outwardly 
increasing ([/' > 0), all energy conditions are violated. If it is introduced at R = a > R max (U is outwardly decreasing, 
V < 0), all energy conditions are satisfied. 

In the regularized geometry, the charge integral (10) involves the contribution due to the charge density on the shell 
and due to the charge density of the ECD distribution outside of the regularizing shell. The resulting (regularized) 
charge density can be written, 

ir° eg = crS(R-a) +pQ{R-a) = s ( R ~ a ) ~ ®{R~a) (26) 

and when integrated over the t = const, slice the expected value m Q 0 is obtained. The Tolman-Whittaker mass 
integral (11) involves the expression 


(2T%-T%) reg 


U' 

4-7T U 3 


S(R — a) + 


U’ 2 

4ttZ7 4 


V 2 t/ \ 
Air U 3 ) 


@(i? — a) 


(27) 


where the first term is due to the energy-momentum contribution of the shell and the second term is due to the 
exterior part of the original ECD distribution. The integral of this quantity over the t = const, slice yields moo, as it 
should be. Thus, the relation (12) holds in the regularized system. 


IV. DISCUSSION AND CONCLUSIONS 

We obtained a spectrum of asymptotically extremal Reissner-Nordstrom (ERN) solutions in the Majumdar- 
Papapetrou (MP) system. The solutions are induced by the assumed distribution of the extremally charged dust 
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(ECD), they have positive ADM mass, the metric profile function U(R) is finite at the harmonic radius R = 0, 
but it can have one or more nodes for R > 0. The outermost node corresponds to the naked pointlike singularity. 
Such solutions are shown to exist only for a discrete set of values of a parameter related to the source strength 
(allocated ECD/mass). Here, for instance, the solutions with zero to three nodes have source strengths: 1.00, 7.29, 
19.6, 39.1. The situation resembles a quantum system with fixed boundary conditions, in which to the normalized set 
of eigenfunctions there corresponds a set of discrete eigenvalues with higher ones corresponding to higher excitations 
(and increasing number of nodes). 

For a distant observer, the ADM mass of the MP spacetime equals the (±) total electrical charge of the system. 
In the same time, the mass of the regular system can be given by the Tolman-Whittaker mass (11), and the total 
charge by the integral of the charge density (10). The spectrum of solutions with nodes in the metric profile function, 
due to the presence of the singularity, does not satisfy the expected relation among the different mass and charge 
parameters (12). Also, since the singularities are naked, the solutions can be taken as physically unacceptable. We 
conclude the discussion of these solutions by showing that by surrounding the singularity with a spherical (5-shell of 
ECD the solutions can be regularized. Interestingly, depending on the position of the (5-shell, the energy density of 
the ECD can have either a positive or a negative sign (see also [16, 17]). For solutions with positive ADM mass and 
the outermost node is at R = Rq, at some R = R m ax > Rq there is a maximum of U(R). If the (5-shcll is placed in the 
region R £ (i?o, -Rmax) the energy density of the shell is negative, while if it is placed in the region R £ (i? max ,oo), 
the energy density is positive. In the special case where the shell is placed exactly at R = R max , the required amount 
of ECD on the shell is zero, but in all these cases the required relation between mass and charge which one finds in 
the regular MP systems is restored. In the similar context, (5-shclls were used to regularize the multiple black hole 
ERN spacetime in [18]. However, (5-shell regularization is not the only possible approach; thick shells [3, 14] could 
nonetheless be applied to achieve the similar results. Alternatively, instead of altering the ECD distribution to replace 
the central singularity with a patch of the regular spacetime, one could, along the lines of Ref. [19], also consider 
using the distributional techniques to construct the distributional energy-momentum tensor that allows the curvature 
singularity to be included in the manifold. 
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